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Abstract 

In this paper, we study the consensus problem for networked dynamic systems 
with arbitrary initial states, and present some structural characterization and di- 
rect construction of consensus functions. For the consensus problem under similar 
transformation, we establish some necessary and sufficient conditions by exploiting 
the structure of consensus functions. Finally, we discuss the consensus problem for 
dynamic systems under switching by using the common Lyapunov function method. 

Keywords: Consensus problem, structural characterization, constructive method, ar- 
bitrary switching, similar transformation, average consensus. 

1 INTRODUCTION 

In recent years, the study of synchronization and coordination of multi-agent systems has 
attracted many researchers. It has broad applications in cooperative control of unmanned 
aerial vehicles, scheduling of automated highway systems, formation control of satellite 
clusters, and distributed optimization of multiple mobile robotic systems ( |T] - pil] ) . 

Consensus problem has a long history. On many occasions, a group of dynamic agents 
in multi- agent /multi-robot systems need to reach an agreement on certain quantities of 
interest. For example, flock of birds tends to synchronize in migration in order to resist 
external aggression and reach their destination. Robots need to arrive at agreement so as 
to accomplish some complicated tasks. Investigation of such problems is of significance 
in theory and in practice. 

Consensus problem was introduced and formally stated by [2|-[I|. in the basic 
definitions were given and average consensus problem was studied for networks with both 
switching topology and time delays. In this paper, we generalize the consensus problem 
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and formulate it in a more general form. |3j tackled this problem mostly by graph theory 
and assumed that the state of each agent is a real scalar. However, in most cases, the 
quantities of each agent are very complex and many aspects should be considered. For 
example, the quantities might be position, velocity, temperature, momentum, voltage, 
mass, energy and so on. Furthermore, these quantities might not be independent. Hence, 
it is natural to extend the domain R of the state of each agent to R m . Therefore, all 
the original definitions for consensus problem should be modified correspondingly. In this 
paper, this kind of consensus problem is studied by using linear algebra theory (j2HI22]) 
as basic tool, and some interesting structural characterizations are established. 

This paper is organized as follows: In section II, we generalize consensus problem 
and establish some necessary and sufficient conditions for a linear dynamic system that 
solves a consensus problem with arbitrary initial state. In section III, we focus on the 
structural characterizations of consensus functions and present a simple and constructive 
method to obtain consensus functions. Furthermore, a necessary and sufficient condition 
for a dynamic system that solves the average consensus problem with arbitrary initial 
state is given. In Section IV, the consensus problem under similar transformation is 
discussed. In Section V, the systems that solve a consensus problem under arbitrary 
switching([o'| |§| ITUj) are characterized. Finally, we summarize our main contribution in 
Section VI. For convenience, some concepts and results in graph theory are given in the 
Appendix. 

2 PRELIMINARIES 

In order to introduce the generalized consensus concept, we consider the following linear 
dynamic system: 
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where G R m , Ay G R mxm , i, j = 1, 2, • • • , n. System can be written in the following 
form 

x = Ax, (2) 

where A = [Ay] and x = [x±, x\, • ■ ■ , x^] T . 

We say Xi and Xj agree if and only if Xi = ^(component-wise). Let n = {1, 2, • • • , n}. 
We say system has reached a consensus if and only if Xj = Xj for all i 7^ j, i, j G n. The 
common value of Xi (i — 1, 2, • • • , n) is called the group decision value. Let \ : R mn — > R m 
be a function of n vectors Xi,X2, ■ ■ ■ ,x n (xi G R m ) and x(0) denote the initial state of 
the system. We say dynamic system solves the X' consensus problem if and only if there 
exists an asymptotically stable equilibrium * T ] T of system (0) satisfying 

x i = x{ x {fy) e R m for all % G n. The function \ is called consensus function. The 
special cases of x( x ) — Ave(a;) = (J^iLi x i)l n i x( x ) = niax(xj), x( x ) = min(xj) are called 
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average- consensus, max- consensus and min- consensus, respectively, due to their broad 



applications in distributed decision making for multi-agent systems. If we have x* = x* 
for all i 7^ j, i,j G n, and x* only relies on initial state x(0), we say that the system solves 
a consensus problem. 

Here, we are interested in the system x = Ax which solves the x- consensus problem 
for some consensus function x an d for any x(0) G H mn . For such systems, there are 
some necessary properties. For example, if x G N(A), then x = 1 <g> b, where N(A) is 
the null space of A, 1 = [1, 1, • • • , 1] T G R n , <S> is the Kronecher product, and b G R' m 
is a constant vector. Furthermore, for any initial state x(0), the solution of the system 
converges asymptotically to some equilibrium. Denote the range (column space) of A by 
R(A). We have the following theorem. 

Theorem 2.1. System (CP solves a consensus problem for any initial state x(0) if and 
only if R{A) = R(A 2 ) and each eigenvalue of A is or has negative real part. Moreover, 
if is an eigenvalue of A, then for any x G N(A), there exists a vector b G R m such that 
x = 1 ® b. 

Proof. Necessity. If R(A) ^ R(A 2 ), then we have rank(A) > rank(v4 2 ). Hence, there 
exists a vector y G H mn , y G N(A 2 ) but y ^ N(A). The solution of system (JTJ) with the 
initial value y is x 



At 

e y 



(I mn + At + ^f- + ■ ■ -)y = y + tAy, where /„ 



is identical 

matrix of order mn. Obviously, x does not converge to any equilibrium when t — > oo, 
which is a contradiction. Therefore R(A) = R(A 2 ). 

Sufficiency. If all eigenvalues of A have negative real parts, then the system is asymp- 
totically stable and all solutions converge to 0, i.e., it solves a consensus problem. If is 
an eigenvalue of A, then there exists an invertible matrix T by R(A) = R(A 2 ), such that 



A = T- 1 
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where J2, • • • , J s are Jordan blocks, and the eigenvalue of Jj has negative real part. So 
the system converges asymptotically to some equilibrium, and since for any x G N(A), 
there exists a vector b G R m such that x = 1 <g) b, the system solves a consensus problem 
for any initial state. □ 

Corollary 2.1. If system (QJj solves a consensus problem for any initial state, then 
dim(N(A)) < m. 

Corollary 2.2. If system solves a consensus problem for any initial state, then 

K mn = N{A)®R{A), 
where © is the operator of direct sum, and 



A(R(A)) = R(A), A(N(A)) = {0}. 



Proof. From Theorem I2.1[ we obtain rank(A) = rank(A 2 ). The remained proof is trivial. 

□ 

In order to investigate system (P) more insightfully, Q can be formulated in the 
following form 
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where C u , D tj G R mxm such that £)J =1 Aj = for all i,j G n. 

Theorem 2.2. System (QJ) solves a consensus problem for any initial state x(0) if and 
only if 

dimN(A) = dimiV(A 2 ) = dimiV([C£, C£, ■ ■ • , C^f) 
and each eigenvalue of A is or has negative real part. 
Proof. We only need to prove that the condition 

dimN(A) = dunN(A 2 ) = dimiV([C 1 T 1 , Cj 2 , • • • , C^J) 

is equivalent to the conditions that R(A) = R(A 2 ) and for any x G N(A), there exists a 
vector b G R m such that x — 1 ® b. 

Obviously, dimiV(A)=dimiV(,4 2 ) is equivalent to R(A) = R(A 2 ). 

(a) Suppose that dimN(A) = r, and for any x G N(A), there exists a vector b G R m 
such that x — 1 ® b. The equation Ax = must have r linearly independent solutions, 
which implies that there are r linearly independent vectors b\, 62, • • • ,b r G R m such that 
A(l £g> bi) = for any i G {1, 2, ■ ■ • , r}. Substituting l®hi into (jlj), we obtain Cubj = 
for any i G n,j G {1,2, ••• ,r}. Therefore dimiVflCft, Cf 2 , • • • ,CjJ T ) > r. But if 
dimiV([Cj[, C22, • • • , Cnn] T ) > r ' then the number of linearly independent solutions of the 
equation [Cj^C^, • •• ,CjJ T a; = is more than r, which implies that dim(iV(A)) > r, 
which contradicts our assumption. Therefore dimiVQC^, C^, • • • , Cj n ] T ) = r. 

(b) If dimiV([Cj[, C^, • ■ ■ , C*^ n ] T ) = r, then there are r linearly independent solutions 
bi, b 2 , ■ ■ ■ ,b r G R m of equation [C^, C^, ■ • • , C^ n ] T x = 0. Thus the equation Ax = has 
r independent solutions 1 <8> b\, 1 <8> 6 2 , • • ■ , 1 <8> 6 r . Since dimiV(A) = r, we obtain that 
for any x G -AT (A), there exists a vector 6 G R m such that x — 1 (g> 6. □ 

3 THE STRUCTURE OF CONSENSUS FUNCTION 

It is important to have clear understanding of the structure of consensus function in 
studying consensus problem. Hence, in this section, we study the consensus function and 
present some characterizations. 



3.1 Consensus Function is a Time- invariant Quantity 



We still consider system (0). If it solves a consensus problem for any initial state, i.e., it 
satisfies the conditions in Theorem 12 .![ then for Wx G R mn , 36 G R m such that lim e 
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1®6 = x*, x* = [J m , 0, • • ■ , 0] (1(g) 6) = [I m , 0, • ■ ■ ,0] lim e Ai x, where I m is identical matrix 

t— >oo 

of order m. Let = l/mj 0, ■ • • ,0] lim ex, then the system solves the x- consensus 
problem. It is easy to see that the consensus function is determined by A. Hence, if 
system (jH) solves a consensus problem for any initial state, it must solve the x-consensus 
problem for some consensus function x- 



If lim e At x = x* is an equilibrium for any x G R m ™, we have Ax* = 0, i.e., A lim e 
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for any x G R r ' 
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0. Since e At A = Ae At , we have lim e At A 



At 



t—>oo 



^At 



Thus A lim e 

t^oo 

A lim e At = 0. Hence dx j^ = [I, 0, • • ■ ,0] lim e At Ax = 0. So the consensus function x( x ) 
is a time-invariant quantity. (Note that lim e At is a constant matrix.) 

Remark 3.1. If system (QJ) solves the x-consensus problem for any initial state, then 
the consensus problem can not be max- or min- consensus. This is obvious by x( x ) = 
[I, 0, • • • , 0] lim e At x. 

t— >oo 



3.2 A Method to Obtain the Consensus Function 

For a given system, the consensus function can be obtained by calculating lim e At . When 

t— >oo 

all eigenvalues of A have negative real parts, it is easy to obtain that xi x ) = 0- However, 
if is an eigenvalue of A, the calculation of lim e At might be very complex. In what 

follows, we will illustrate that, for some special cases, we can find a simple method to 
obtain the consensus function. 
Consider the following system: 
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denoted by A = Ax, which satisfies the conditions in Theorem l2.1l and rank(A 2 ) 
{n — l)m. 

By Theorem 12 .2[ it is easy to show that 
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where B {j G R mxm , B,; G R mxmn for any i,j G n. 

Since system (jSJ) solves the x-consensus problem for any x(0) G R mn , and let 
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respectively, we get B u = B 21 = ■ ■ ■ = B nl , B 12 = B 22 = • ■ ■ = B n2 , ■ ■ ■ , B ln = B 2n = 
■ ■ ■ = B nn , i.e., Bi = B 2 = ■ ■ ■ = B n . We denote B t by E = (E%, E 2 , • • • , E n ), where 
Ei G R mxm for any % G n, so 



B 



Since lim e At A = A lim e At = 0( mn)x(mn ), we get BA = AB = 0( mn)x ( mn) . Therefore 
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Since xi x ) is an invariant quantity, we have 



lim x(x(t)) = x( 
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This implies 



Ex(0) 

(Ex + E 2 + --- + E n )Ex(0) = Ex{0) 



for all x(0) G K mn . 

By the theory of Jordan canonical form, we learn that rank(i?) = m, i.e., raxik(E) = m, 
so {£x(0)|Vx(0) G R m "} = R m . Hence 



Ei + E 2 + ■ • - + E n = I„ 



(7) 



Because rank(A) = (n — l)m, there exist m linearly independent vectors £i, £2, • ' ' > £n 
in R mn such that ([A = for all i G {1, 2, • ■ ■ , m}. Let 



[^i,Z 2 



where G R mxm . Then there exists an invertible matrix T G R mxm such that = TZ. 



By (JZJ), we have 
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T(Z 1 + Z 2 + --- + Z n ) = I n 



So Z\ + Z 2 + ■ ■ ■ + Z n is invertible and T = (Zi + Z 2 + ■ ■ ■ + Z n ) 1 . Therefore 

x (x) =Ex = (Z 1 + Z 2 + --- + z n y x zx. 
By the discussion above, we get the following procedure to get the consensus function: 
Method 1. 1. Choose arbitrarily m linearly independent vectors £1, £2, ■ • • , Cm G N(A T ); 

2. Let F = (ft, ft, • • • , U) T = {Fi, F 2 ,'--, F m ), where F t G R mxm , i = 1, 2, • • • , m; 

3. Let T = Fx + F 2 + ■ ■ • + F m , then T is invertible; 



4- x( x ) = T 1 Fx, lim e 



At 



T- X F 
T~ X F 

T^F 



3.3 An Example 

In what follows, we present an example to show the effectiveness of Method 1. 
Example 3.1. 



A 



-10 1 
10-10 
10-1 
10-1 



By Theorem 12. 1\ it is easy to verify that system x = Ax solves a consensus problem 
for m — 2. 
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On the other hand, calculating lim e At directly, we get 

t^oo 
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Hence the consensus function obtained by Method 1 is correct. 



3.4 Average Consensus Problem 

Average consensus problem has been discussed in [3] and j3j, and the authors of them 
presented some necessary and sufficient conditions. In this section, we also consider the 
problem but from another viewpoint. 

Based on the discussion in Subsection A, we set average consensus function x( x ) = 
F m x(mn)% = \( x i + x 2 + • ■ • + x n ) such that FA = 0, where x G R mn , X{ G R m for any 
% G n. 

Let 
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respectively, we get 
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Therefore dim(iV(y4 r )) > m. By Corollary 12.11 we have the following lemma. 

Lemma 3.1. If system (Op solves the average consensus problem, then rank(A) = (n — 
l)m. 

Theorem 3.1. System (QJ) solves the average consensus problem if and only if rank(A 2 ) =ran 
m(n — 1), [/,/,■•• , I] A = and A[I, I, ■ ■ ■ , I] T = 0, where I G R mxm is identical matrix, 
and is an eigenvalue of A, and all the other eigenvalues have negative real parts. 



Proof. The proof is obvious. We omit the details. 



□ 



3.5 The Case of m = 1 



In this subsection, we study the consensus problem in the case m = 1. 
Consider the system 



x = Ax, 
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and Xj, a^- G R, Vi, j G 

Theorem 3.2. System (0j solves the \- consensus problem for any initial state if and only 
if 

1 ) each eigenvalue of A is or /ias negative real part; 

2) if zs an eigenvalue of A, then Al = 0, and rank(74 2 )=rank(yl) = n — 1. 



So we get a general method to derive the consensus function. 
Method 2. Choose arbitrarily y G N(A T ), y ^ 0, and let y = [2/1,2/2 



, y n ] T , then the 



consensus function is 



where x = [x\, ■ ■ ■ , x r , 



En 
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We have the following corollary. 

Corollary 3.1. System (0) solves the average consensus problem if and only if all the 
nonzero eigenvalues of A have negative real parts, rank(A 2 ) =rank(A) = n — 1, 1 T A = 0, 
and Al = 0. ([S], Theorem 5) 

Remark 3.2. Naturally, if A is a Laplacian matrix of some graph, then it satisfies the 
conditions in Corollaru \S.l\ except 1 T A = 0. 

Remark 3.3. we can view the consensus function 

En 
i=\ Ui x i 
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as a weighted average consensus function. 



4 CONSENSUS PROBLEM UNDER SIMILAR TRANS- 
FORMATION 



Since consensus function is determined by A, the study on the structure of A becomes an 
important issue. 



Let 
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where M is a nonsingular real matrix, the eigenvalues of M have negative real parts and 
r < m. Hence, by Corollary 12.21 if system (0) solves the x-c° nsensus problem for any 
initial state, and dim(N(A)) = r, then there exists a nonsingular real matrix T such that 
A = T~ l JT. We will study the structure of T in this section. 

From the discussion in Section II, system (0) solves the x- consensus problem for any 
initial state if and only if lim e At exists and the equilibriums have the form of 1 ® b. In 

t^oo 

the following, we will show that, for some systems, if lim e At exists, we may find similar 

n^oo 

transformation T such that T~ X AT solves the x-consensus problem. 
Consider the following system 
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Theorem 4.1. System (0) with r = m solves the x-consensus problem for any initial 
state if and only if 
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where Tij G R mxm for any i,j G n, satisfies Ty = Y^j=iTij is invertible, Y^j=iTij = 
for i = 2, ■ • ■ , n, and 
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is invertible. 



Proof. Sufficiency. It suffices to prove that all equilibriums of T l JT have the form: 
1(g) b. 
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, b m G R m be linearly independent vectors, and let Ci 
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1 ® b m , which are also linearly independent, then T JTQ = 0, % 



1, 2, • ■ ■ , m. The sufficiency is proved. 

Necessity. The equation T~ x JTx = must have m linearly independent solutions: 



Cl = 1® 61, C2 = 1®&2 



Cn 



where 6j G R m , 
Notice that 



, m. 



T~ l JTx = <^> JTx = 0^Tx 



c 








for some c G R m . So there exist 


Cl, C2, • • • 


, c m G R m such that 








c 2 




Cm 

















TCi = 








j ■ ■ ■ j T(^ m 











n 

u 




n 


and 


























1 










( £?*)&* = 


0. 


i = 2,3, ••• ,n, 








1 









for k — 1, 2, • • • , m. Since &i, 62, ■ • • , &m are linearly independent, we have 



^^. = 0,i = 2,3, 



n. 



Since T is invertible, we obtain that Ty and To are invertible. 

For the average consensus problem, we have the following theorem. 



□ 



Theorem 4.2. System solves the average consensus problem for any initial state if 
and only ifr — m and 
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where Y^=i Tij = 0, % = 2, 3, • • • , n, and Tu and To all are invertible. 

Proof. Sufficiency. Let 
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[/,/,-•• ,I]T- 1 JT = 0. 

By Theorem 13.14 T~ X JT solves the average consensus problem. 

Necessity. If system © solves the average consensus problem, we have 



T~ l lim e Jt T 
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We obtain that Tn = Tn = ■ ■ ■ = T\ n are invertible. 

By Lemma 13.14 we get r = m, and by Theorem 14. 1| we have YTj=i Tij — 0, i 
2, 3, • • • , n and T D is invertible. 

For the case of r < m, T is very complex, but we still have the following theorem. 

Theorem 4.3. If 
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where Tjj and Tp defined as in Theorem J^.l are invertible, and X]j=i = 0, 
2, 3, • • • , n, then system solves the x-consensus problem for any r < m. 



Proof. For arbitrary r linearly independent vectors b\ , b 2 , 



b r G R r , we define 
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mxl , 
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for z = 1,2, • • • , r. Then (1,(2,- 
1,2,-.. ,r. 

Therefore, system solves the x-consensus problem. □ 

Remark 4.1. //lim^oo e A * exists and dim(N(A)) = r < m, then we can find an invertible 
matrix T such that system x = T~ l ATx solves the x-consensus problem. 

Moreover, if r = m, we can find T such that system x = T~ x ATx solves the average 
consensus problem. 

Example 4.1. For any initial state, the system in Example \3. 1\ solves a consensus prob- 
lem, but not the average consensus problem. We will provide the procedure to find invert- 
ible T such that system x = T~ l ATx solves the average consensus problem. 
First, we can choose an invertible matrix T\, 
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C = T 2 7 1 BT 2 = 

By Theorem 13. 1[ x = Cx solves the average consensus problem for any initial state. 



Let 
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then x = T 1 ATx solves the average consensus problem for any initial state. 



5 CONSENSUS PROBLEM UNDER ARBITRARY 
SWITCHING 



In this section, we investigate the consensus problem of system 

A(t)x. 
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The study on the consensus problem of (|TU|) is difficult. Here, we only consider some 
special cases. We view (J10J) as a switched system and A(t) is a constant matrix in each 
switching interval. 

We consider the following system, each subsystem of which is the same as 



x(t) = A s( t)x(t), 



(11) 



where s(t) : R + — > n is the switching signal. 

Generally speaking, not all switched systems solve a consensus problem for any initial 
state. But some special switched systems can solve a consensus problem. 

We assume that, for Vs G n, 

1) x = A s x solves a consensus problem; 
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where A S ij is a symmetric and positive definite matrix (denoted by A S ij > 0) for all 
3) A s a = 



x 



A s ij -i 

4) every subsystem has the same consensus function x( x ) = F m 
where F { e R mxm , F, > for all i G n, and FA S = 0; 

5) FiA S ij = A S ijFi for all i, j G n, which implies FiA S ij > 0. 
Before presenting Theorem 15.11 we first prove the following lemma. 
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Lemma 5.1. Let 
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be a symmetric matrix, where Lij G R mxm ; Lij > and YTj=i j^i Lij = La for alli,j G n. 

If the eigenvalues of L are arranged in an increasing order Ai < A2 < • • • < \ mn , then 
we have 

Xi = X 2 = ■ ■ ■ = X m — 0, A m+ i > 0. 
Then we call L a block laplacian matrix. 

Proof. For any x G R mn , we have 

n 

X 1—jX — ^^^^ X J-j ^X ^ ^^^^ X J-j^jX j 

8=1 i^j 



— ^ (x^ LijXi ~\~ Xj LijXj Xj LijXj Xj LijXij 



2 ^ j i x i x j) Lij(xi Xj) 
9 ^ ~~ ^i) Lij(xi — Xj). 



Therefore Ai = A 2 



0, A m+ i > 0. 



□ 



Theorem 5.1. If switched system Ml)) satisfies the conditions (l)-(5), then, under arbi- 
trary switching, the solution of the system globally asymptotically converges to l®x{ x {^))> 
i.e., the switched system solves the x-consensus problem. 

Proof. Since x( x ) is an invariant quantity for every subsystem, xi x ) is a lso an invariant 
quantity under switching. For any solution x(t), let x(t) = 1 £g> x{ x ) + 5(t). We refer to 5 
as the (group) disagreement vector. Then 



Fx = F(l <g> x(x)) + FS = Fx + FS. 



Hence, we have 



F5 = 0. 

For a given s G n, since x = A s x, we have 



(13) 



dt 



+ 5 = A,(l ® x(x)) + A,5, 



which implies 
Let 



5 = A J. 



e 



Fi 



then 6 > 0. 

Let V{8) = 5 T Q5, then 

^ = 25 T Q5 = 5 T (BA S + A^Q T )5. 

Let —L = QA S + Aj6 T , then L is a block laplacian matrix by assumption. We can 
easily get that N(L) = R(l <g> J) and N(L) n iV(F) = #(1 ® I) fl iV(F) = {0}. We divide 
the linear space R mn into the direct sum of and its orthogonal complement space 

R(l ® J)- 1 , then we have 



R mn = ^(1 0/) 0j R(l 0/) 



Correspondingly, 5 = 5\ + 5 2 , Si G R(l <S> I), £2 G R(l <S> -0 -1 . Let P be the orthogonal 
projector from R' mn onto R^l^I) 1 - such that 62 = PS. Since F5 = 0, we have P5 = 62 7^ 

if <y ^ 0. 

Hence <f = -5 T L5 = -(<? x + 5 2 ) T £(5i + 6 2 ) = -5jP T LP5 2 < -\ m+1 5 2 r P T P5 2 < 0, 
where A m+ i > is the (m + l)th smallest eigenvalue of L. This shows that V(S(t)) is 
a valid common Lyapunov function for the group- disagreement, i.e. , under arbitrary 
switching, the switched system solves the x-consensus problem. □ 

Remark 5.1. The assumptions (1-5) seem rather strict, but this kind of system really 
exists extensively. For example, it is easy to show that the system 

x = -Lx, (14) 

where L is block laplacian matrix, solves the average consensus problem and satisfies the 
assumptions (1-5). 

6 CONCLUSIONS 

For linear dynamic systems, consensus problem has been discussed from a new viewpoint. 
The structure of the consensus functions has been characterized. An example has been 
presented to illustrate the effectiveness of our results. Some necessary and sufficient condi- 
tions for consensus problem under similar transformation have also been obtained. Finally, 
we characterize a class of dynamic switched systems that solve a consensus problem under 
arbitrary switching. 

7 APPENDIX: GRAPH THEORY PRELIMINAR- 
IES 

In this section, we briefly summarize some basic concepts and results in graph theory that 
are useful in dealing with the consensus problem. More comprehensive discussions can be 
found in [23] . 

A undirected graph Q consists of a vertex set V = {ni,n 2 , • • • ,n m } and an edge set 
£ = {(jii,rij) : Ui,nj G V}, where an edge is an unordered pair of distinct vertices of V. 
If ni,rij G V, and (ni,rij) G 8, then we say that and nj are adjacent or neighbors. An 
oriented graph is a graph together with a particular orientation, where the orientation of 
a graph Q is the assignment of a direction to each edge, so edge (rij, rij) is an directed edge 
(arc) from 77^ to rij. The incidence matrix B of an oriented graph Q is the {0, ±l}-matrix 
with rows and columns indexed by the vertices and edges of Q, respectively, such that 
the r/'-entry is equal to 1 if edge j is ending on vertex rii, -1 if edge j is beginning with 
vertex rij, and otherwise. Define the Laplacian matrix of Q as L(Q) = BB T . L{Q) is 
always positive semi-definite. Moreover, for a connected graph, L(Q) has a single zero 
eigenvalue, and the associated right eigenvector is l m . 
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